KRL-MAP-309 



Cosmological Signatures of Interacting Neutrinos 

Nicole F. Bell, 1 ' 2 -E Elena Pierpaoli, 2 'Q and Kris Sigurdson 2 - 3 '!] 

Kellogg Radiation Laboratory, California Institute of Technology, Pasadena, CA 91125 
2 Theoretical Astrophysics, California Institute of Technology, Pasadena, CA 91125 
!1 School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540 

(Dated: March 29, 2006) 

We investigate signatures of neutrino scattering in the Cosmic Microwave Background (CMB) 
and matter power spectra, and the extent to which present cosmological data can distinguish be- 
tween a free streaming or tightly coupled fluid of neutrinos. If neutrinos have strong non-standard 
interactions, for example, through the coupling of neutrinos to a light boson, they may be kept in 
equilibrium until late times. We show how the power spectra for these models differ from more 
conventional neutrino scenarios, and use CMB and large scale structure data to constrain these 
models. CMB polarization data improves the constraints on the number of massless neutrinos, 
while the Lyman-a power spectrum improves the limits on the neutrino mass. Neutrino mass limits 
depend strongly on whether some or all of the neutrino species interact and annihilate. The present 
data can accommodate a number of tightly-coupled relativistic degrees of freedom, and none of the 
interacting-neutrino scenarios considered are ruled out by current data — although considerations 
regarding the age of the Universe disfavor a model with three annihilating neutrinos with very large 
neutrino masses. 



I. INTRODUCTION 

We are in a remarkable era when cosmological data 
is both precise and abundant. Measurements of tem- 
perature and polarization fluctuations in the cosmic mi- 
crowave background (CMB) radiation (e.g. 0, HE ||), 
the matter power spectrum via galaxy surveys [Jj, la and 
the Lyman-a forest 0, @ , the current and past expan- 
sion rate of the Universe via the Hubble Key Project 
and observations of Type la supernovae |lfj| respectively, 
and the abundance of light elements predicted by big 
bang nucleosynthesis (BBN) all show striking con- 
sistency within the standard ACDM cosmological model. 
Robust bounds can be placed on the forms of matter 
and energy that constitute the Universe, and cosmol- 
ogy is now a powerful particle physics laboratory that 
constrains the properties of both the new dark-matter 
particles demanded by cosmological data (e.g. 01) an( i 
the familiar particles of the standard model of particle 
physics. In many cases these constraints would be diffi- 
cult or impossible to obtain in any other way. 

Neutrino physics is an excellent example. Cosmologi- 
cal techniques for probing neutrino properties rely upon 
detecting indirect signatures of the relic neutrino back- 
ground, which complements the significant experimental 
and theoretical effort that is underway to understand the 
surprising physics of the neutrino sector. A particular fo- 
cus of this effort has been to set limits on the mass and 
number density of relic neutrinos. In this work we in- 
stead focus on the cosmological signatures of neutrino 
interactions. 
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It is well known that the CMB can be used to con- 
strain the number of light relativistic degrees of freedom 
(in addition to the photon) - conventionally parametrized 
as the effective number of neutrino species, Nf. A rel- 
ative increase in the radiation density delays the epoch 
of matter-radiation equality, and leads to an enhanced 
Integrated Sachs- Wolfe (ISW) effect. Present CMB lim- 
its are 1.6 < Nf < 7.1 [ll fl3|. Big Bang nucle- 
osynthesis (BBN) also constrains Nf, as additional ra- 
diation increases the expansion rate and alters the ex- 
pected primordial abundance of helium. The current 
BBN bound is (up to various flavor dependent subtleties) 
Nf < 3.3-4 [13. 

The standard cosmological model predicts Nf « 
3.04, consisting of the three known neutrino species, 
plus a small correction that accounts for the neutrino 
heating from electron-positron annihilation and finite- 
temperature QED effects (e.g. |l5|). With only the three 
standard neutrino species, BBN constraints, combined 
with neutrino mixing, no longer permit the possibility of 
a significantly enhanced Nf due to large chemical poten- 
tials |16(. However, it is important to bear in mind that 
Nf may include not only neutrinos, but any light par- 
ticles that are thermally populated in the early universe. 
In models containing sterile neutrinos or other light rel- 
ativistic degrees of freedom, the BBN limits may be sub- 
stantially modified, and Nf > 3 is still possible ^tJ- (I n 
addition, Nf < 3 can be obtained if the reheating tem- 
perature following inflation is low 0|.) Moreover, the 
relativistic energy density may evolve between the time 
of BBN and CMB decoupling, so Nf\mN_a,nd A^ ff |cMB 
need not be the same quantity (e.g. jl9l l20j| ) . 

In addition to the total relativistic energy density, cos- 
mology can also be used to probe interactions in the neu- 
trino sector. Neutrinos with only standard model cou- 
plings interact via the weak force, and decouple when 
T ~ 1 MeV. At later times, they free-stream, and inter- 



act only via gravity. In this paper we investigate a class 
of models which feature extra, non-standard, neutrino in- 
teractions. In these models, neutrinos interact strongly 
with a new scalar boson, which is brought into thermal 
equilibrium though its coupling to the neutrinos. Rather 
than free-streaming, the neutrinos form a tightly coupled 
fluid with the new scalar. 

These models generically have non-standard values for 
N° s , but perhaps more interestingly, the absence of neu- 
trino free-streaming leaves a distinctive signature in the 
CMB. If the neutrinos are part of a tightly coupled 
fluid, they are fully characterized by density and ve- 
locity perturbations, and anisotropic stress is negligible. 
In |2lL 122 ] it was shown that the current Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) CMB measurements 
already have some sensitivity to this effect. This is sig- 
nificant because in addition to being able to infer the 
presence of relativistic degrees of freedom, we may now 
also be able to say something about the interactions of 
the particles which make up that relativistic energy den- 
sity. 

In this paper we address the question: how much rel- 
ativistic energy density is there, and what fraction of it 
must consist of weakly interacting particles ? We answer 
this question in general, and also in the context of specific 
models. 



II. INTERACTION MODEL 

Although the results of our analysis are valid in a wider 
context than the interaction model we now describe, we 
examine in this section a simple physical model of non- 
standard neutrino interactions for illustrative purposes. 

We consider the coupling of neutrinos to each other 
with bosons, through tree level scalar or pseudo-scalar 
couplings of the form 



C-vtj) = hijViVj(j> + gijUi^Ujcf), 



(1) 



where the boson (j> is taken to be light or massless 1 . Such 
couplings arise in Majoron-like models, viable examples 
of which have been discussed in Ref . [2_^| . Recently, these 
models have been investigated in the context of late-time 
phase transitions, whereby the neutrinos acquire their 
masses via a symmetry breaking phase transition at a 
low scale, which occurs late in the history of the universe 
[l^. [2{| . In order to be as model independent as possi- 
ble, we assume the new couplings are fixed independently 
of the neutrino mass. We also make no distinction be- 
tween g or h type couplings, nor between neutrinos and 
antineutrinos. 

Existing bounds on these new couplings are extremely 
weak. For example, the solar neutrino [26] and meson 



FIG. 1: The interactions that keep the neutrinos and the 
scalar coupled. If the scalar is heavier than m„, the process 
v <-> v<j) is replaced by cj> <-> vv. 



decay [23 limits are \g\ < 10~ 2 . Neutrinoless double 
beta decay sets a limit g ee < 10~ 4 28], but does not 
constrain other elements of the coupling matrix g a p. 
Supernova constraints exclude a narrow (and model- 
dependent) range of couplings around g ~ 10~ 5 pfjj . 
Even couplings which are much smaller than these limits 
can have significant cosmological consequences. 

For a massless (j> boson, scalar couplings could medi- 
ate long-ran ge forces with possible cosmological conse- 
quences |30L l31|. while pseudo-scalar couplings mediate 
spin-dependent long-range forces, which have no net ef- 
fect on an unpolarized medium 2 . However, if the <j> boson 
has even a tiny mass Hq <C <C 1 eV the interaction is 
short ranged and insignificant over cosmological distance 
scales. 

The <p boson can be brought into thermal equilibrium 
through its coupling to the neutrinos, and the v — <p sys- 
tem may stay in thermal contact until late times. The 
processes involved, shown in Figure ^ are v 4> v 4>i 
vv 4>(f>, vv vv, and either v +-> vcf> or vv <-> 0, 
depending on whether the scalar mass, m^, is smaller or 
larger than the neutrino mass, m u 3 . For sufficiently large 
couplings, the v— <j> system will remain in thermal contact 
until the temperature drops below m„ or m<f,. At this 
point the heavier of the two particles will annihilate or 
decay. 

The possibility of altering the relativistic energy den- 
sity through neutrino decay has been considered in |33| 4 . 



2 For pseudo-scalar couplings, two-boson ex change can mediate 
extremely weak spin-independent forces |32| . 



1 Couplings of neutrinos to new heavy bosons are tighty con- 
strained l23l. 



3 We set all three neutrino species to a common mass m„, with 
niu ^» \J ^^goi i y <^ m atm • When this approximation does not 
hold, the effects of neutrino mass are negligible in present cos- 
mological data; 

4 See also, Ref. l'S4l . which studies the case of a scalar boson de- 
caying into neutrinos, thus distorting the usual thermal neutrino 
distribution. Related scenarios, in which hot dark matter is pro- 
duced by the decay of heavier particles, are examined in Ref. l35l . 
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while the cosmological effects of neutrino annihilation or 
self- interaction were examined in [3||[32]]. In particular, 
[3fil | considered the introduction of self-interactions as a 
mechanism to eliminate neutrino free-streaming, thus ob- 
taining neutrinos which behave as cold dark matter, de- 
spite their light mass. However, in these scenarios the 
neutrinos were taken to be heavy enough (m v ~ lOeV) 
to contribute all of the dark matter, which is no longer a 
viable possibility. As we will show below, the combined 
use of neutrino mass limits and cosmological observations 
now allow much more sensitive constraints to be placed 
on these types of models. 

The ultimate effect of the interaction shown in Eq. (JJJ 
can be split into two distinct cases, depending on whether 
the scalar is lighter or heaver than the neutrinos. We now 
discuss these cases in turn. 



A. Light Scalars, < m„ 

When TO0 < m„, the neutrinos can annihilate to 
scalars via the process vv — > (fxj) when the temperature 
drops below m„. Complete annihilation occurs provided 
that g > 10 -5 , although smaller couplings would suffice 
to keep the v—<j> system in thermal equilibrium until late 
times via decay/inverse decay processes. Obtaining par- 
tial annihilation (rather than cither negligible annihila- 
tion or complete annihilation) would require fine-tuning 
of g. 

If all three neutrino species completely annihilate, they 
will make no contribution to the dark matter density in 
the universe today, as we would be left with a "neutrino- 
less universe" at late times. This eliminates cosmological 
constraints on neutrino mass poj . However, a nonzero <fi 
mass would make a small contribution to the dark mat- 
ter density. Note that < m v implies the neutrinos 
are unstable, and may decay via i>i — > Vj<j>. For the same 
range of couplings, g > 10~ 5 , this could lead to neutrino 
decay over astronomical distances, which is testable in 
future neutrino telescope experiments |38j . 



B. Light Neutrinos, m$ > m v 

Alternatively, if ra^ > m v , the scalar will eventually 
convert to neutrinos, so that cosmological neutrino mass 
bounds do apply. For example, the decay/inverse-decay 
process <f> <-> w, will keep the v— 4> system tightly coupled 
throughout the CMB era, provided g > 10 -12 |19|. Once 
the inverse process becomes kinematically inaccessible at 
T ~ 777a, the 4> bosons will decay to neutrinos. 



III. TIGHTLY COUPLED PERTURBATION 
EQUATIONS 

For most parameters of interest, the neutrino mean free 
path is much smaller than the scales of interest during 



the epoch of CMB decoupling. In this limit, the neutri- 
nos and bosons form a single, tightly coupled fluid. Per- 
turbations in this fluid evolve differently, compared to 
the usual collisionless neutrino background. As they are 
coupled via gravity, neutrino perturbations influence the 
evolution of photon perturbations and thus neutrino per- 
turbations can leave a distinctive signature in the CMB. 
In this section we discuss the relevant properties of the 
v-(j) fluid for interesting limits of the model considered in 
Section ITTland in a more general context. 

In the standard scenario, free-streaming damps the 
neutrino density perturbations, and introduces a source 
of anisotropic stress, e.g. see [3!l|4C|. In comparison, a 
tightly-coupled fluid has only density and velocity pertur- 
bations, with the shear stresses and all higher moments 
in the Boltzmann hierarchy absent (at least to linear or- 
der). Defining the density and velocity perturbations as 
8 = dp/ p and 9 — ikv respectively, the equations describ- 
ing the evolution of the tightly-coupled v-<j> fluid are 




-(! + <") 



.« (1 _ 3w )0--^- 
a K ' 1 +( 



-k 2 S, 



(2) 
(3) 



in the synchronous gauge, where an overdot is a deriva- 
tive with respect to conformal time. Here, c 2 = SP/8p is 
the adiabatic sound speed and u> — P/p is the equation 
of state with p = p v + p$ and P = P v + P$. In the 
limit where both v and <j) ar e relativistic ui = c 2 = 1/3. 
However, for non-zero masses, u> and c 2 temporarily de- 
crease from 1/3 and deviate from each other, during the 
period when v or <fi starts to become non-relativistic and 
annihilate/decay. We may define an effective number of 
standard model neutrinos N™ 1 contributed by the v-<p 
fluid with the relation 



7tt 2 
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(4) 



where = (4/H) 1 / 3 ^ is the canonical cosmic neu- 
trino background temperature in the standard model. 
The value of A^ nt after the heavier species annihilates 
into the lighter will be greater than the value before an- 
nihilation. Details of the calculation of w, c 2 , and N™* 
can be found in Appendix lAl 

The differences from standard cosmology in these mod- 
els will therefore typically be a combination of (i) elimina- 
tion of neutrino free-streaming, (ii) nonstandard equation 
of state and sound speed evolution, and (iii) additional 
(and perhaps an evolving amount of) relativistic energy 
density. We now consider two phenomenological models 
(limiting cases of the interaction model of Section[Hj that 
exhibit these differences. 
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A. Model A: Freestreaming vs. Interacting 
(Massless Particles) 

In this section we consider a scenario in which some 
fraction of the neutrinos act as a tightly-coupled fluid, 
but where N° s is constant with time. We parameterize 
this class of models by N^ M and -/V™*, the number of 
standard model (free-streaming) and interacting degrees 
of freedom respectively, where 



Nf = iV, SM 



(5) 



A scenario in which N° s is constant in time corresponds, 
for example, to the limit m v — > and — > 0. However, 
even for finite masses, this limit is a good approximation 
as long as N% s does not evolve significantly during the 
times of interest, as would be the case if any annihila- 
tion/decay takes place well before or well after the CMB 
decoupling epoch. In particular, this is a reasonable ap- 
proximation for the 77i0 > m v models considered in [19j, 
where only modest evolution of occurs. 

Note that this description is much more general than 
any particular limit of a neutrino-scalar interaction 
model. It encompasses any scenario in which some frac- 
tion of the energy density is in free-streaming relativistic 
particles, and another fraction is in a tightly-coupled rela- 
tivistic fluid, as long as N° s does not evolve significantly. 

Fig. El shows effects on the CMB angular power spec- 
trum and the matter power spectrum of varying the num- 
ber of interacting neutrinos, while iV® and all other cos- 
mological parameters are held fixed. This allows us to 
explore the effects of suppressed neutrino free streaming 
alone. In the standard case (A™* = 0) when a perturba- 
tion of a given scale enters the horizon, power is trans- 
fered from the neutrino density modes to higher moments 
of the neutrino distribution as the neutrinos free-stream 
out of gravitational potentials. This effect does not occur 
for interacting neutrinos which instead contribute to the 
gravitational potential and thus enhance the monopole 
perturbation of the photon distribution and increase the 
amplitude of the CMB temperature power spectrum for a 
fixed amplitude of primordial perturbations. Indeed, we 
see that beyond the first peak the effect of increasing A" 1 * 
in the CMB can be roughly approximated as a constant 
offset in the amplitude of the spectrum. A small shift in 
the location of the peaks is also a feature . In the up- 
per panel of Fig.[3]we show the contributions to the CMB 
power spectrum arising from the monopole terms (the 
physical temperature perturbation at the last-scattering 
surface corrected for its gravitational redshift), the veloc- 
ity terms (perturbations due to the Doppler shift), and 
the ISW terms (perturbations due to the evolving gravi- 
tational potentials) of the source function (see, for exam- 
ple, Refs. |4j, |42| for a pedagogical review of the source 
function). The lower panel of Fig. shows the difference 
in CMB power spectra between models with AT™* = 3 and 
A™* — keeping Nf = 3 fixed (the extreme models of 
Fig- El- While the change in the monopole dominates the 
total difference in the power spectrum at all I, the ISW- 
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FIG. 2: CMB and matter power spectra as a function of the 
fraction of interacting neutrinos, with Nf = N^ M + Nl nt = 3. 
The power spectra are normalized (to an arbitrary value) at 
large scale. 



monopole cross term contributes significantly to the dif- 
ference near the first peak. Note that this ISW-monopole 
contribution is nonzero solely because gravitational po- 
tentials evolve differently in a model with and without 
free-streaming neutrinos — the background evolution is 
identical in the two cases. 

As we are in the limit of massless neutrinos, the effect 
on the matter power spectrum of suppressing free stream- 
ing is very minor. Free-streaming has a significant effect 
on the matter power spectrum when m„ is finite, such 
that neutrinos contribute some fraction of the dark mat- 
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FIG. 3: Upper: The CMB power spectra due to the 
monopole, velocity, and ISW terms of the source function 
for the standard Nl nt = case (Fo r a p edagogical description 
of these terms see, e.g., Refs. |4ll li^l. Lower: The contri- 
bution of each source term to AC;, the difference between a 
model with Nl nt = 3 and Nl nt = with Nf = 3 held fixed. 
The ISW-monopole cross term (nonzero due to the differing 
evolution of the gravitational potentials) contributes signifi- 
cantly to the difference in the first peak despite the identical 
background evolution of the two models. 



ter density today. In that case, the free-streaming of the 
neutrino hot dark matter component damps the growth 
of structure while the neutrinos are still relativistic. By 
comparison, the effect of modifications to the neutrino 
perturbations of massless neutrinos is very small. 



FIG. 4: The CMB and matter power spectra as a function of 
the number of standard model neutrinos, with Nl nt = 0. The 
power spectra are normalized (to an arbitrary value) at large 
scale. 



For comparison, we show in Fig. 0] the effects on the 
CMB angular power spectrum and the matter power 
spectrum of varying the total relativistic energy density, 
while N™ 1 and all other cosmological parameters 
are held fixed. We see that increasing enhances the 
first peak. This is the result of a larger ISW effect, due to 
the delay in matter-radiation equivalence. There is also 
a large shift in the positions of the subsequent peaks, 
which occurs due to the change in the conformal time of 
last scattering. See Ref. for a discussion of these ef- 
fects. In Fig. we also see that the enhanced radiation 
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FIG. 5: The evolution of the effective number of interacting 
neutrinos with m„ = 0.1 eV (dashed curves), m„ = 1 eV 
(dotted curves) and m„ = 10 eV (solid curves) for 3 in- 
teracting neutrinos (top/red), 2 interacting neutrinos (mid- 
dle/green), and 1 interacting neutrino (bottom/blue). Notice 
that the Nl nt initially includes an extra 4/7 to account for the 
new scalar degree of freedom and that the effective number of 
interacting neutrinos after annihilation is greater than prior 
to annihilation. 



density suppresses the matter power spectrum. Again, 
this is because matter-radiation equivalence, and hence 
the growth of structure, is delayed. 

The extent to which the various effects discussed above 
can be compensated with a change in other cosmological 
parameters will be discussed in section Ivl 



B. Model B: Neutrino Annihilation to Scalars 

We now turn to the slightly more complicated scenario 
with nifp < m„, where the interacting neutrinos can an- 
nihilate. We again assume the limit — > 0, but allow a 
nonzero m u . We specialize to the case of three neutrino 
species, and will consider the possibility that either one, 
two, or all three neutrino species interact strongly with 
the scalar </>, which we shall denote by models Bl, B2, and 
B3 respectively. For couplings constants g > 10~ 5 , the 
neutrino species which are coupled to the scalar will an- 
nihilate when T v ~ m v . If all three neutrino species an- 
nihilate (model B3) this leaves a "neutrinoless universe" . 
Cosmological neutrino mass bounds are altered in these 
scenarios, because the neutrino species which annihilate 
will not make a contribution to the dark mater density 
today (i.e., they will not contribute to fl u ) pfjj . 

A distinctive feature of this scenario is that the neu- 
trino annihilation will heat the scalars, causing N* to 
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FIG. 6: The evolution of w, the equation of state (solid 
curves) and (? a , the sound speed (dotted curves) as a function 
of the scale factor, a, for m v = 1 eV for 3 interacting neutrinos 
(bottom/red), 2 interacting neutrinos (middle/green), and 1 
interacting neutrino (top/blue). Both deviate from the value 
of 1/3 for a relativistic fluid, during the epoch of annihilation. 



evolve as the annihilation proceeds. For simplicity, we 
shall assume the scalar boson is brought into thermal 
equilibrium before the neutrinos thermally decouple from 
the electrons and positrons (at T ~ 1 MeV). The scalar 
will then initially contribute an amount 8N V = 4/7 to the 
relativistic energy density, so that N° s ~ 3.57 5 . This 
value will increase as the annihilation proceeds. For ex- 
ample, if all three neutrinos annihilate, the final relativis- 
tic energy density is equivalent to N v ~ 6.6 [23] (see also 
Appendix El . For realistic neutrino masses, m v ~ — 2 
eV, this annihilation occurs close to the time at which 
the CMB photons last scatter, so that the evolution of 
takes place during the CMB decoupling era. The 
evolution of N° s is shown in Fig. The evolution of the 
equation of state, w, and sound speed, c~, during this 
annihilation epoch is shown in Figure. [5] 

The present laboratory limit on neutrino mass is m v < 
2.2 eV, set by tritium beta decay experiments 0. Given 
the tiny mass squared differences measured by solar and 
atmospheric neutrino oscillation experiments (i5m 2 ol ~ 
7 x 10~ 5 eV 2 and <5m 2 tm ~ 2 x 10~ 3 eV 2 0), the tri- 
tium bound applies to all three neutrino mass eigenstates. 
We shall assume that the three neutrino eigenstates have 
degenerate masses, which is a good approximation for 



5 If the scalar were not populated until sometime after the neutri- 
nos thermally decouple from the e+e - plasma, N 1 ^ would not 
be altered (as energy density would simply be shifted from one 
relativistic species to another.) 
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FIG. 7: The CMB and matter power spectra, for varying 
values of m v , for model Bl (one interacting neutrino and two 
standard neutrinos). The power spectra are normalized (to 
an arbitrary value) at large scale. 



FIG. 8: The CMB and matter power spectra, for varying 
values of m„, for model B3 (three interacting/annihilating 
neutrinos). The power spectra are normalized (to an arbitrary 
value) at large scale. 



m v <L 0-1 e V- (We use m„ to denote the value of a 
single neutrino mass throughout, so that the quantity 
Ysm v = Sm^.) In the analysis in section we consider 
neutrino masses in the allowed range m v = — 2.2 eV. 

For simplicity, we also assume that the scalar has suf- 
ficiently strong self-interactions that it continues to be- 
have as a tightly-coupled fluid once the neutrino annihi- 
lation is complete. For early neutrino annihilation (large 
masses) this will result in the largest deviations of the 
CMB spectra with respect to the standard scenario. For 
late annihilation (small masses) the late-time behavior of 
the scalar is irrelevant. 



Fig. shows effects on the CMB angular power spec- 
trum and matter power spectrum of varying m„ in the 
model with one interacting neutrino (BI), while Fig. 
shows the same effects for the model with three inter- 
acting neutrinos (B3). Here we see a combination of the 
effects of free-streaming suppression (compare FigEJ) and 
larger N° (compare Fig.^J. In the CMB spectra, we see 
the enhanced overall amplitude which is characteristic of 
the tightly coupled neutrinos, together with an enhanced 
first peak and shifted subsequent peaks that result from 
increased relativistic energy density. As expected, the ef- 
fects on the CMB spectra are more pronounced in model 
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B3, in which all three neutrinos species annihilate. More- 
over, the effect increases as m u increases, because larger 
mass corresponds to earlier annihilation, so that the ex- 
tra relativistic energy density is present earlier. 

Conversely, the effects on the matter power spectra 
are more pronounced in model Bl in which only one 
neutrino species annihilates. This is expected, because 
the remaining neutrinos contribute to the dark matter 
density today. This remaining neutrino hot dark mat- 
ter component causes the usual suppression of the power 
spectrum that is used to constrain neutrino mass. How- 
ever, in the "neutrinoless" model (B3), this suppression 
of the power spectrum is absent, because no neutrino 
dark matter component remains. A smaller power spec- 
trum suppression does remain, as shown in Fig. [S] This 
small effect is a result of enhanced N° s , which slightly 
delays matter-radiation equivalence and hence delays the 
growth of structure. Again, this effect is more significant 
for larger m v , because the extra radiation is present ear- 
lier. We thus see that the "neutrinoless universe" model 
(B3) can accommodate large neutrino mass , wh ile having 
little effect on the matter power spectrum [20| . 



V. RESULTS 

We now discuss how the different models fit the data. 
We first note that for all of the models, parameters sets 
can be found which provide a good global fit to the data, 
and so none of the models can be ruled out. For example, 
in Tabled we report the x 2 values of the best fit points 
in our Markov chains for the case of the Cosmo data set. 
The table shows that the best fit parameter set for all 
the models considered indeed provides a good fit to the 
data. We plot spectra for the best fit models in Fig. [5] 
Apart from the high value of N„ in model A (see the 
discussion at the end of the the next section) , the best- 
fit parameters are within commonly adopted parameter 
ranges in the standard ACDM model. 



Model 


x 2 


DOF 


Reduced \ 2 


(DOF/2)- 1 / 2 


A 


1515.00 


1478 


1.025 


0.037 


Bl 


1520.17 


1479 


1.028 


0.037 


B2 


1523.99 


1479 


1.030 


0.037 


B3 


1526.77 


1479 


1.032 


0.037 



TABLE I: The \ 2 values and de grees of freedom (DOF) , using 
the Cosmo data set. 



IV. DATA 

We computed the CMB and large scale structure power 
spectra, and performed the parameter estimation with a 
modified version of the publicly available Markov Chain 
Monte Carlo package, COSMOMC [U[46|. 

Constraints on the models were evaluated using CMB 
data from the WMAP % ACBAR and CBI Q ex- 
periments, together with the galaxy power spectra mea- 
sured by the 2dF Galaxy Redshift Survey [j| and the 
Sloan Digital Sky Survey (SDSS) @|. In addition, we 
used the measurement of the Hubble parameter made by 
the Hubble Space Telescope (HST) Key Project [§]. This 
set of data constitutes our basic comparison set (hereafter 
Cosmo). For model B3, we imposed a cut-off of 2.2 eV 
for the neutrino mass, consistent with the tritium beta 
decay limit ■ 

In addition, we will investigate how the introduction of 
the Lyman-a constraints Q and the CMB polarization as 
measured by CBI Q constrain these models. These two 
variants will be named CosmoLya and CosmoCBIpol re- 
spectively. We implement the Lyman-a constraints de- 
rived in Ref. Q and we implement them in a similar way 
as in Ref. 8], with minor modifications that were sug- 
gested by the authors. 

Note that SDSS and 2dF data impose a constraint only 
on the shape of the power spectrum, P(k), and not on 
the normalization, as we make no assumptions about the 
bias. Lyman-a data, however, imposes constraints on 
both normalization and shape at small scales. 



On the other hand, while reasonable parameters can 
be found that yield a good global fit to the data for each 
of model A, Bl, B2, and B3, we will show that mod- 
els with fewer interacting neutrinos are preferred in a 
Bayesian sense after marginalizing over all other cosmo- 
logical parameters. The interpretation of these results is 
discussed further below. We shall now discuss each model 
individually, focusing on the constraints on the neutrino 
properties. 

A. Model A: N* M vs. 7V; nt 

Fig- El displays the curves for the marginalized likeli- 
hood of the parameters N„ and iV" 1 * . Notice that the 
CMB alone allows a larger number of interacting neu- 
trinos, while preferring a relatively low number of stan- 
dard neutrinos. The addition of the matter power spectra 
yields the effect of reducing the maximum value of iV^ nt , 
and also shifts the peak of the likelihood for N^ M . Re- 
ferring to Fig. |21 we notice that keeping the total number 
of neutrinos fixed while increasing the number of inter- 
acting neutrinos has a large impact on the CMB power 
spectrum, and a very minor one on the matter power 
spectrum (the two models have equal total energy den- 
sity, but different neutrino fluctuation evolution). How- 
ever, the introduction of the matter power spectrum con- 
straint in the likelihood analysis imposes limits on the 
epoch of equivalence, and on the total spectral index. 
This, in turn, leaves less freedom for accommodating 
non-standard neutrino fluctuation evolution which af- 
fects the CMB power spectrum. 
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FIG. 9: The CMB temperature (TT, upper left), polarization (EE, upper right) and cross (TE, lower left) power spectra, and 
the matter power spectrum (lower right), for the best-fit models. In each case models can be found that fit all data well. 



As anticipated in |14|. the addition of the Lyman-a 
constraint improves the limits on the standard model 
neutrinos. This is because by probing the power of 
the fluctuations on very small scales, the Lyman-a data 
places a strong constraint on the spectral index n @, 
which is degenerate with N^ M and the reionization r in 
the radiation power spectrum. 

It is interesting to note that adding the CBI polariza- 
tion data also improves the constraints on Nf . This 
is because an increased number of neutrinos significantly 
shifts the peaks of the Ci, due to the delay in matter- 
radiation equivalence. The CBI polarization data have 
been shown to be able to determine the phase of the os- 
cillations with high precision Q| , despite the size of the 
crrorbars. This is one example of how such information 



can be used in constraining parameters. 

Table ITT| summarizes the best fit values and the 95 % 
confidence levels for the marginalized likelihoods for N^ M 
and Nl nt , while Fig. 1111 shows the degeneracy between 
N„ and N™ 1 in the case of the Cosmo dataset. A higher- 
value for A^ nt may be compensated by a lower value for 
N„ M . However, once the matter power spectrum is used 
in the analysis (and thus the redshift of equivalence is 
constrained) the degeneracy is mild. 

It is clear from Figs.llUI andllllthat models with fewer 
than three interacting neutrinos are favored in a Bayesian 
sense. This does not contradict our previous statement 
that cosmological parameter values can be found with 
iV^ nt ~ 3 which provide a good global fit to the data 
considered (a reduced \ 2 statistically consistent with 
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FIG. 10: The marginalized likelihood curves for N' v nt and N° M 
in Model A. The solid curves corresponds to the Cosmo data 
set (CMB+2dF+SdDS+HST) and the long-dashed curve to 
the CMB data alone. The dotted and short-dashed curves 
correspond to the addition of the Lyman-a and CBI polariza- 
tion data respectively. 







N SM 


Nl nt 


Nl nt 


Data set 


best fit 


upper limit 


best fit 


upper limit 


Cosmo 


5.0 


7.8 


0.015 


2.9 


CosmoCBIpol 


4.3 


6.9 


0.1 


2.9 


CosmoLya 


4.9 


6.8 


0.02 


3.0 



TABLE II: The best fit values and 95 % C.L. allowed ranges 
for 7V^ M and in Model A. The confidence limits are ob- 
tained from marginalized curves. 



1). It is worth elaborating on the implications and in- 
terpretation of this Bayesian limit and why no contra- 
diction exists. Firstly, the likelihoods shown are after 
marginalizing over all other parameters assuming flat pri- 
ors in N„, Nl nt , and all other cosmological parame- 
ters. As such, the likelihood in Fig. ^]is not straight- 
forward measure of the goodness-of-fit of models with 
given values of (iV^ M ,A^ nt ) but also of the density of 
good- fitting models with the values (N^ M ,Nl nt ) in the 
hyperplane of parameter space defined by the remain- 
ing cosmological parameters — it technically quantifies 
the cumulative relative likelihood of two populations of 
models. We can certainly conclude, for instance, from 
Fig. ^] that (with the aforementioned flat priors) the 
set of data we have used prefers the population of mod- 
els with (N^ M ,Nl nt ) ~ (3,0) relative to the population 
of models with (N^ M , N™ 1 ) ~ (0,3) at more than 2-cr. 




0.5 1 1.5 2 2.5 

J\Jtot 



FIG. 11: Upper: The contour level for the parameters N„ 
and Nl nt in Model A. The two solid lines are one and two 
a contours, shades are from the mean likelihood, using the 
Cosmo data set. Lower: The relative likelihood along the 
line Nt nt + = 3. 



However, we can not conclude from Fig. llll that all mod- 
els with (N^ M ,Nl nt ) ~ (0, 3) are ruled out at more than 
2-cr by current data. 6 This conclusion can only be made 



The notion of "more than 2-<r" depends on the context. In refer- 
ence to the marginalized likelihood in a two dimensional param- 
eter plane it refers to a given point in this parameter plane being 
outside the 95% confidence region. In reference to a particular 
model with a given set of cosmological parameters it refers to 
that model having a reduced x 2 more than two standard devia- 
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if no models with (N^ M ,N'l nt ) ~ (0,3) can be shown to 
be consistent with the data, which as we have discussed, 
is not the case. Similarly, while we can conclude that 
the population of models with (A^ M , j\P nt ) ~ (3,0) are 
disfavored compared to the population of models with 
{N^ M , AT™*) ~ (5, 0) we can not conclude that all models 
with 5 free-streaming neutrinos provide better global fit 
to the data than all models with 3 free-streaming neu- 
trinos. This is just the nature of Bayesian inference — 
it attempts to quantify the relative likelihood of models 
within a given paradigm, but does not make an abso- 
lute judgment on the viability of the paradigm itself (or 
in this case a particular subset of the models within the 
paradigm). To make a judgment on the viability of a 
given paradigm, a better approach is to ask the ques- 
tion "Are there any models within this paradigm which 
are consistent with the data?". If the answer to this 
question is "yes" then that paradigm is still viable. In 
this case, another way of summarizing the state of affairs 
is that even if priors stipulating (A^ M ,i\P nt ) ~ (3,0) 
or (A^ M ,A^ nt ) ~ (0,3) are imposed to restrict our 
paradigm to either the standard model of particle physics 
or an alternative model with interacting neutrinos, then 
models that are globally consistent with the data can still 
be found (albeit fewer of them in the later case). 

Finally, let us comment on parameter degeneracy. Be- 
cause the number of neutrinos contribute to setting the 
redshift of equivalence, both N^ M and A^ nt are de- 
generate with U, m and Hq. It has been pointed out 
|l3l Il4| that a higher number of massless neutrinos can 
be compensated by a higher Hq. We find that even the 
marginalized likelihood for Hq in this model is around 80 
km s _1 Mpc -1 , which is only one a away from the HST 
quoted best fit. In this respect, our results are discrepant 
with those in |'2l| , where much larger values of Hq were 
obtained. 



B. Model B: Annihilating neutrinos 

In Fig. E| we plot the marginalized likelihood for the 
neutrino mass in model Bl (one interacting neutrino plus 
2 standard model neutrinos) using the Cosmo and Cosmo 
+ Lyman-a data sets (CosmoLya). The curves for B2 
are similar to Bl. In Fig ll31 we show the corresponding 
likelihood for the case where all three neutrinos annihi- 
late (model B3). The best fit values and the 95 % C.L. 
of the marginalized likelihoods are reported in table Hill 
for the model in which either one, two or three neutri- 
nos interact and annihilate. For models Bl and B2 the 
best fits are at m„ ~ eV. Note that our approximation 
does not allow us to explore m v ~ eV, as the assump- 
tion that all three neutrino masses are equal breaks down 
when m v < y/ <5m^ tm ~ 0.05 eV. For model B3 the best 



tions (of the \ 2 distribution) away from 1. 




FIG. 12: Marginalized likelihood for the neutrino mass in 
Model Bl. The solid curve and dashed curves correspond to 
the Cosmo and CosmoLya data sets respectively. 




FIG. 13: Marginalized likelihood for the neutrino mass in 
the "neutrinoless universe" model, B3. The solid and dashed 
curves correspond to the Cosmo and CosmoLya data sets 
respectively. Note that we imposed a 2.2 eV cut-off for m v . 

fit is non-zero, and is discussed further below. 

For the cases Bl and B2 where only one or two neutrino 
species annihilate, the remaining neutrinos contribute to 
the dark matter density today. As expected, Lyman- 
a data significantly tightens the neutrino mass limit in 
these models, because it sets a constraint on the over- 
all normalization and the shape of the power spectrum. 
This reduces the 95 % C.L. upper limit by quite a lot, 
bringing it to 0.24 eV for one interactive neutrino (model 
Bl) and 0.31 eV for 2 interactive neutrinos (model B2). 
The addition of the CBI polarization data does not im- 
prove the constraints on neutrinos mass with respect to 
those obtained with the Cosmo data. This is because 
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Bl 


Bl 


B2 


B2 


B3 


B3 


Data set 


best fit 


limit 


best fit 


limit 


best fit 


limit 


Cosmo 


< 0.1 


1.5 


< 0.1 


0.42 


2.2 


2.2 


CosmoLya 


< 0.1 


0.24 


< 0.1 


0.31 


2.2 


2.2 



TABLE III: The limits on the neutrino mass, in the various 
models. All values are in eV. For each model, the upper limit 
refers to the 95 % C.L. limit of the marginalized likelihood. 
For m v < 0.1 eV, the approximation that all three neutrino 
masses are degenerate is not satisfied. Note that for model 
B3, we imposed a 2.2 eV cut-off, as implied by the tritium 
beta decay neutrino mass limit. 



such low neutrino masses don't have a significant impact 
on the CMB power spectrum (see Fig. 0). As for degen- 
eracies, the neutrino mass parameter is degenerate with 
il m and Ha and Hq, all of which affect the redshift of 
equivalence and therefore the amplitude of the matter 
power spectrum at small scales. As a consequence, m v is 
also degenerate with the amplitude of the matter power 
spectrum as- The mean of the marginalized Hq likeli- 
hood plot is 68 km s" 1 Mpc -1 and 71 km s" 1 Mpc -1 for 
models Bl and B2, which are perfectly normal values. 

The case B3 is different, because here all three neutri- 
nos species annihilate, leaving no neutrino contribution 
to the dark matter today. The effects on the power spec- 
trum are thus milder (see Fig. [HJ) , and hence the limits on 
the neutrino mass are expected to be significantly weaker. 
In particular, the more massive the neutrino is, the ear- 
lier it annihilates (see Fig.|5jl. If the neutrinos annihilate 
early enough that all cosmologically interesting scales are 
outside the horizon, then models corresponding to differ- 
ent masses only differ by a slightly different expansion 
history very early on. Therefore, the data cannot distin- 
guish between different high values of m„, as the corre- 
sponding power spectra are very similar (see Fig. |H1 and 
Fig. El) 

However, for neutrino masses vn v ~ 1 eV, the anni- 
hilation takes place very close to the time of CMB de- 
coupling. It is during the annihilation period that the 
greatest deviation of the sound speed and the equation of 
state occur (see Fig. [HJ) Hence for m v ~ 1 eV, the CMB 
spectrum is affected by the modified values of w and c 2 s , 
in addition, of course, to the lack of free-streaming. This 
tends to disfavor m u ~ 1 eV, with respect to both larger 
and smaller values of vn v , as shown in Fig. 1131 

In addition, the neutrino mass in model B3 shows sig- 
nificant degeneracies with the Age and Hq, as changes in 
these parameters compensate for the effect of increased 
radiation density on the epoch of equivalence. The neu- 
trino mass is also degenerate with tti,h 2 and as- We have 
plotted these degeneracies in Fig. El To help explore 
the degeneracies, we have determined the marginalized 
likelihood contours with and without imposing the 2.2 
eV cut-off in m u . We can see from Fig. El that although 
a slightly better fit is obtained for larger m„ (5-10 eV), 
this would imply an unacceptably low Age |47j . an un- 




FIG. 14: The confidence level contours for the parameters 
flbh 2 , Age, Ho, and as, vs. m v , for model B3. The solid lines 
are the one and two sigma contours, and shades are from the 
mean likelihood, using the Cosmo data set. The panels on the 
left have a 2.2 eV cut-off imposed on m„, while the panels on 
the right have a 15 eV cut-off. 



acceptably high Hq, and also a value of erg > 1, which 
is disfavored by recent cluster number counts |48j and 
lensing analysis |49j |. (The baryon abundance however, 
more closely matches the BBN determination, fl^h 2 = 
0.022 ±0.002 |Jj|> f° r larger values of m v .) Imposing the 
2.2 eV cut-off for m v (as required to be consistent with 
the tritium beta decay neutrino mass limit) brings these 
parameters back to more reasonable values. However, 
even in the 0-2.2 eV region, the values of these param- 
eters seem somewhat discrepant. For example, the best 
fit point (at m v = 2.2 eV) has values of f^/i 2 , Age and 
Hq which all deviate from the central values preferred by 
BBN, globular cluster, and HST measurements, respec- 
tively, by about one sigma. 
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C. Discussion 

The results above should be compared with the ones 
obtained by other authors in the standard case of three 
massive, non-interacting neutrinos. The WMAP team 
found m v < 0.23 eV at 95 % CL [U, using CMB data 
and the 2dF power spectrum, plus an assumption on 
the power spectrum bias which drives most of the re- 
sult [5l], |52] • The addition of quite tight constraints from 
the SDSS power spectrum leads to an upper limit on the 
sum of the three masses of VJ m„ = 0.42 eV at 95 % 
CL [5i| . The recent results from the Boomerang team 
are m v < 1 eV at 95 % CL from CMB data alone, and 
m v < 0.4 eV when the matter power spectrum data are 
considered, with no bias assumptions |54j. Perhaps the 
most directly comparable value is the one obtained by 
[55}, ^rriv — 0.33 ± 0.27 eV, with the same Lyman-a 
sample and comparable assumptions about the parame- 
ters as in our analysis. These results are broadly similar 
to what we find for the cases Bl and B2, in which one or 
two standard model neutrinos remain. As a general re- 
sult, we can conclude that the neutrino mass is expected 
to be below about 0.5 eV in those cases, and is not sig- 
nificantly affected by the presence of the tightly-coupled 
component of the relativistic energy density. However, 
the "neutrinoless" model (B3) is different, and leads to 
much weaker m„ constraints, albeit at the expense of 
somewhat discrepant values for ilbh 2 , Age and Hq. 

Finally, we comment on the implications of these re- 
sults for neutrino decay. One consequence of strong v — <p 
interactions is that neutrinos may be unstable. For the 
range of couplings considered here, g > 10 ~ 5 , the neu- 
trinos could decay over astronomical distances, which is 
testable in neutrino telescope experiments |3^] . For small 
masses, the cosmological consequences of these couplings 
are represented by Model A, while for larger masses they 
are represented by either Model B2 or B3 (depending on 
whether one or two neutrino species are unstable.) 

Ref. [5(j claims stringent limits on the couplings g^, 
and thus on neutrino lifetimes, based upon the fact 
that Refs. 0, find some evidence for neutrino free- 
streaming in the CMB. However, 0,H2 analyzed only a 
scenario in which all three neutrinos behave in the same 
way, and thus certainly cannot be used to claim that all 
three neutrino species must be free-streaming. 7 

Given the results here, we conclude that current cos- 
mological data does not impose limits on either the indi- 
vidual couplings, gij, or upon neutrino lifetimes. 



Neutrino decay requires only some subset of the couplings to be 
large (e.g. 921 allows for the decay ^2 — > z^i + </>) implying that 
only a subset of the neutrinos, Ui, need be interacting. Without 
demonstrating that all the couplings g^ must simultaneously be 
small, it is not possible to set stringent cosmological limits on 
any particular individual element of gij. 



VI. CONCLUSIONS 

We have investigated the distinction between a free- 
streaming or interacting fluid of neutrinos, and explored 
the constraints on these neutrinos imposed by cosmology. 
Although some of our results are applicable in a wider 
context, we have used an example involving additional 
couplings between the neutrinos and a light boson. Our 
main conclusion is that models with interacting neutrinos 
remain viable, contrary to the claim in Ref. pH [5|| . 

As a general result, we find that both CMB polariza- 
tion data and Lyman-a data help to constrain the num- 
ber of standard model neutrino species. The inclusion 
of Lyman-a data also helps improve the constraints on 
neutrino mass. With the data considered, we have found 
upper limits on the neutrino mass, but no detection of a 
nonzero neutrino mass. 

Two parameterizations of the interacting-neutrino 
models have been examined: 

(1) In the first, we allowed for an arbitrary number 
of free-streaming (standard model) neutrinos, A^ M , and 
tightly-coupled (interacting) neutrinos, N™ 1 . Within the 
context of the neutrino-scalar model, this corresponds 
to the limit m v — > and — > 0. We have found 
that within the (N^ M ,Nl nt ) plane, the data favors free- 
streaming neutrinos over tightly-coupled neutrinos in a 
Bayesian analysis. However, we find even if the prior 
jySM is imposed (so that all neutrinos are interact- 
ing) models can be found that are a good global fit to 
the data considered, without resorting to extreme values 
for the cosmological parameters. In this respect, models 
with interacting neutrinos remain viable. We empha- 
size that the constraints obtained here are very general, 
as N„ and JV™* can parametrize any free-streaming 
or tightly coupled relativistic degrees of freedom, which 
need not consist of neutrinos. 

(2) Our second parametrization consists of models with 
non-zero m v . We fixed the total number of neutrinos to 
three, and allowed cither one, two, or all three, to interact 
with a massless boson. In this scenario, the interacting 
neutrinos annihilate to bosons when T ~ m„, thus re- 
moving them from the plasma. In the case where either 
one or two neutrino species interact/annihilate, we find 
the upper limits on neutrino mass are broadly similar 
to those for the standard scenario. However, if all three 
neutrino species annihilate to leave a "neutrinoless uni- 
verse" , the neutrino mass limits are significantly weaker. 
In this case, values of m v comparable to the tritium beta 
decay limit of 2.2 eV are permitted, although a low Age 
and a high Hq tend to somewhat disfavor the scenario. 
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is the Riemann zeta function of three, and the arrows in- 
dicate the high-temperature limit x — * 0. Notice that the 
standard result for massless particles Pi — pi/3 is recov- 
ered in this limit. These results can be straightforwardly 
derived by expanding the distribution function as a ge- 
ometric series of Boltzmann factors. This form is useful 
for tabulating the density and pressure as the K a can be 
rapidly evaluated and the sum converges quickly and can 
be truncated to the desired accuracy — keeping only the 
leading term is the Maxwell-Boltzmann approximation 
to the distribution function. 

Specializing to the fermionic case (£ = 1) we write the 
energy density and pressure of a single fermionic degree 
of freedom as 



APPENDIX A: NEUTRINO-SCALAR FLUID 
PROPERTIES 

Here we summarize a few basic relationships for Fermi- 
Dirac and Bose-Einstein statistics and derive from them 
the quantities w, c s , and iV^ nt used in Eqs. ©-0) to 
evolve the perturbation dynamics and shown in Fig. 
and Fig. El We note that these quantities may also be 
found by solving the continuity equations for energy and 
entropy density numerically, but we believe the analytical 
forms we present here are physically illustrative and may 
be useful in other contexts. 

The number density rii, energy density pi, and pressure 
Pi of a fermionic (£ = 1) or bosonic (£ = —1) species i of 
mass rrii with gi internal degrees of freedom in thermal 
equilibrium can be written as 
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where Xp/p ~~* 1 f° r x ~ * but are < 1 for finite x. We 
also find it useful to define 
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where similarly Xdp/dP * 1 f° r x ~ > 0. 

Let us now consider, as we do in this paper, a thermal- 
ized fluid at temperature T v $ consisting of N massive 
neutrinos coupled to a single massless scalar degree of 
freedom cf>. 

The total energy density and pressure of this fluid are 



- Pu = + 2N~ Xp (x„)T^ (A8) 



and 



P = P <P + P V = ^T* + 2N~ XP (x v )T^ (A9) 



where E = ^/p 2 + mf, x = m,i/T, K a are modified 
Bessel functions of the second kind, £(3) = 1.2020569 



where the factor of 2 in the neutrino terms accounts for 
anti-neutrinos and x v = m^/T^. 



15 



We can now write the equation of state w 
terms of Xp and Xp as 



P/p in fluid is constant. The constancy of S implies that 



4/7 + jVxpQzy) 
4/7+Nxp(x v ) 



Similarly, we can write the soundspeed c 2 
terms of Xdp and XdP as 



4/7 + N Xd p(x v ) 
4/7 + N X dp(xv) 



(A10) 
dP/dp in 

(All) 



The form of Eq. (|A10|) and Eq. IjAlljl explains the 
behavior shown in Fig. El For T u( f, ^> m v we have x ~ * 1 
and both w and c 2 approach 1/3. For T v $ <C m v we 
have x ~ * and again w and c 2 approach 1/3. It is only 
during the annihilation of the neutrinos (T„0 ~ m u ) that 
the values of w and c 2 deviate from those for a relativistic 
fluid. Larger values for iV result in larger deviations. 

We now derive an expression for the temperature of the 
neutrino-scalar fluid T v< j, as a function of time (measured 
with T 7 ). In the standard case the neutrino temperature 
just falls as T$ M = (4/H) 1 / 3 !^ cx a. It is convenient to 
measure T„0 relative to this standard case and define the 
ratio lZ U( p = T„,f,/T^ M . For the times of interest the weak 
interactions have decoupled already and the comoving 
entropy density S = a 3 (p + P) jT v ^ of the neutrino-scalar 



A/7 + N 



4/7 + N[(3/4) Xp (x v ) + (l/4) XP (av)] 



1/3 



(A12) 



which, recalling that x v — m v / (R u< pT^ M ) , is a transcen- 
dental equation which implicitly determines lZ vc t> as a 
function of T™. 

Now if we write the energy density in terms of an ef- 
fective number of standard model neutrinos 



(A13) 



we find, comparing with Eq. 



that 
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where the last line holds for T u( f, <C m v (after annihilation 
is complete). In accordance with Fig. [3] we find that 
NT -> (2.20,4.25,6.58) for N = (1,2,3) respectively. 
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